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Abstract. LS I +61° 303 is one of the most observed Be/X-ray binary systems because, peculiarly, it has periodical radio 
and X-ray emission together with strong, variable gamma-ray emission. This source remains, however, quite enigmatic. Some 
properties of this system can be explained assuming that the unseen companion is a non-accreting young pulsar with a relativistic 
wind strongly interacting with the wind of the Be star. On the contrary, other properties of LS I +61° 303 fit a model where the 
companion is accreting even with two events of super-critical accretion along the orbit. The very recent discovery of a radio 
jet extending ca. 200 AU at both sides of a central core has definitely proved the occurrence of accretion-ejection processes 
in this system. Therefore it is of great interest to combine this result with previous observations at other wavelengths within 
the framework of the two-peak accretion (ejection) model. Concerning the first ejection, we show that the observed gamma- 
rays variations might be periodic with outbursts confined around the periastron passage (i.e. where the first accretion-rate peak 
occurs and high-energy emission but no radio emission is predicted). Concerning the second ejection, with radio bursts, we 
point out that it can be also traced in the X-ray data, both in episodes of hardening of the X-ray emission and in a transition 
from soft- to hard-states at the onsets of radiobursts. In fact, both hardening and transitions between spectral states are related to 
the dramatic change in the structure of the accretion disk preceeding the ejection. Finally, we explore the nature of the accretor 
and we conclude that on the basis of the present optical data a black hole cannot be ruled out. 



Key words, stars: individual: LS I +6F303, 2CG 135+01 - X-rays: binaries 
vations - gamma-rays: theory 



radio continuum: stars - gamma-rays: obser- 



1. Introduction 

The X-ray binary stellar system LS I +61°303 has always ex- 
cited particular interest because of its two peculiarities: On the 
one hand being the probable counterpart of the variable gamma 
ray source 2CG 135+01 (Gregory & Taylor [T978I Tavani et 
al. Qp98 ) and on the other hand being a periodic radio source 
(Taylor & Gregory [13521 [T3g4l 

The binary system is composed of a compact object in 
an eccentric orbit around a rapidly rotating B0-B0. 5 main- 
sequence star, undergoing mass loss through an equatorial disk 
(Hutchings & Crampton ll981> . The orbital period is assumed 
to be of 26.496 days determined on a radio data base of over 
20 years (Hutchings & Crampton ll98U Gregory et al. 1999 
Gregorv l2002l . 

High radio outbursts always peak near phase 0.6, where 
<3>=0 has been set at Julian Date 2443366.775 (Taylor & 
Gregory 1982 Paredes et al. 1990). However, lower in- 
tensity radio outbursts occur within the broad distribution 
<t> (fl( /, o =0.45-0.95, due to variations of the equatorial disk in the 
Be star (Paredes et al. 1 19901 Gregory et al l 19991 Zamanov & 
Marti llUUni Gregory et al. l^UOH Gregory & Neish l^UOlt . From 
optical and near-infrared observations the periastron passage 
is estimated to be in the range ^ perjastron -0.2-Q.5 (Hutchings 
& Crampton [19811 Marti & Paredes [19951 . One of the fun- 
damental questions concerning the periodic radio outbursts of 



LS I +61° 303 therefore has been: Why are the radio outbursts 
shifted with respect to the periastron passage? 

The X-ray emission is also periodic with period estimates 
: 26.7 : 



of P= 26.7 + 0.2 by Paredes et al. ( 1997 1 and P=26.42+ 0.05 by 
Leahy ( 2001 ), clearly in agreement with the radio period. Quite 
surprising, however, is that the X-ray outbursts are offset from 
the radio ones. In fact, the two available simultaneous X-ray 
and radio observations by Taylor et al.( 1996 1 and by Harrison 
et al. (|2000) during an orbital period show that the X-ray emis- 
sion peaks in the phase interval <I>x-ray=0.43-0.47 (recalculated 
by Gregory 2002 1 while the radio outburst is offset to that by 
several days. 

Taylor et al. J 19921 and Marti & Paredes J1995l > have mod- 
elled the properties of this system in terms of an accretion rate 
M oc S^-, (where p W i„d is the density of the Be star wind 
and v re i is the relative speed between the accretor and wind) 
which develops two peaks for high eccentricities: the high- 
est peak corresponds to the periastron passage (highest den- 
sity), while the second peak occurs when the drop in the rela- 
tive velocity v re i compensates (because of the inverse cube de- 
pendence) the decrease in density. From typical parameters of 
LS I +61° 303 derived from near infrared data, Marti & Paredes 
(1995) have shown that both peaks are above the Eddington 
limit and therefore one expects that matter is ejected perpen- 
dicular to the plane of the accretion disk. Near the periastron 
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the short distance to the Be star enhances inverse Compton 
losses: X-ray or/and gamma-ray outbursts are expected but no 
radio bursts. At the second accretion peak, the compact object 
is much farther away from the Be star, so that the electrons can 
propagate out of the orbital plane: an expanding double radio 
source should be observed with a radio interferometer (Taylor 
& Gregory [T984l Taylor et al. [T992l Massi et al. [T993l Marti & 
Paredes [T^95l Gregory et al. [T999l 

The main problem concerning this model is that the lu- 
minosity of LS I +61° 303 in the X-ray range is only Lx — 
10 35 erg s _I (Maraschi & Treves lT98U . That is three orders of 
magnitude lower than the Eddington limit, even for a neutron 
star, and in addition the bulk of the energy output seems to be 
shifted from X-ray to y-ray wavelengths with L y st 10 37 erg s" 1 
(Hartman et al. 119991 . The difficulty of interpreting these re- 
sults in the context of a super-accretion model has led to an 
alternative young pulsar model where a population of relativis- 
tic electrons are produced at the shock boundary between the 
relativistic wind of the young pulsar and the wind of the Be star 
( Maraschi & Treves fTWTl Tavani [T995l Harrison et al. 120001 
Hall et al. 2003)- As a matter of fact such a model fits the time- 
variable high-energy emission observed near periastron from 
the Be/pulsar system PSR B 1259-63 (Tavani & Aaron [T997t . 

However, the recent discovery of a radio emitting jet ex- 
tending ca. 200 AU at both sides of a central core (Massi et 
al. 2004 1 has shown the occurrence of accretion-ejection pro- 
cesses in LS I +61°303. Therefore, this source seems very sim- 
ilar to the microquasar LS 5039 (Paredes et al. 2000 1 also sub- 
luminous in the X-ray range (even more than LS I +61°303) 
and also having L y > Lx- The quite stable gamma-ray emis- 
sion in that case is explained due to upscattered stellar photons 
via inverse Compton from the relativistic electrons of the per- 
sistent jet. If that is true for LS I +61°303 (Taylor et al. 1996 
Massi et al. 2001 , Kaufman Bernado et al. 20021, do the peri- 
odic outbursts of this source imply periodic gamma-ray bursts 
? And in this case at which <b y - ra y ? Moreover, recent devel- 
opments in the theory of the accretion-ejection processes show 
that magneto-rotational instabilities are able to accelerate and 
to collimate a part of the disk material into a double jet only af- 
ter thermal instabilities in the accretion disk have inflated and 
transformed it into a geometrically thick disk (Meier 2001}. As 
shown by Belloni et al. ( 1997 1 structural changes of the disk are 
associated with changes in the X-ray spectral states. Is it possi- 
ble to discern a change of state in available X-ray data ? And in 
this case: are changes of state and radio outburst related to each 
other as expected in an accretion-ejection process ? The aim of 
this paper is to try to answer these questions. Section 2 analy- 
ses and discusses gamma-ray data while Section 3 deals with 
X-ray and optical results. The conclusions are given in Section 
4. 

2. EGRET data analysis 

The EGRET gamma-ray data in Fig. [2 extend over a period 
of 1300 days. During that long interval there are three data 
sets "a,b," and "c" lasting for 10, 7 and 22 days, respectively. 
In order to search for eventual periodicitites we used the Phase 
Dispersion Minimization (PDM) method which is very efficient 



on irregularly spaced data (Stellingwerf 1978 1. We used the 
PDM algorithm of the UK Starlink software package, PERIOD 
( http ://w w w. starlink.rl . ac.uk/ 1 . 

As shown in Fig.|2]the PDM results in a dominant feature 
at P=27.4+7.2. The folding of the EGRET data with this period 
results in a clear clustering of the gamma-ray emission (Fig.|3)- 

The statistical significance of the period is calculated in 
PERIOD following the method of Fisher randomization as out- 
lined in Nemec & Nemec d!985> . The advantage of using a 
Monte-Carlo- or randomization-test is, that it is distribution- 
free and that it is not constrained by any specific noise mod- 
els (Poisson, Gaussian etc.). The fundamental assumption is: If 
there is no periodic signal in the time series data, then the mea- 
sured values (gamma-ray flux in our case) are independent of 
their observation times and are likely to have occurred in any 
other order. One thousand randomized time-series are formed 
and the periodograms calculated. The proportion of permuta- 
tions that give a peak power higher (for the PDM: a trough 
lower) than that of the original time series would then provide 
an estimate of p, the probability that no periodic component is 
present in the data. A derived period is defined as significant 
for p<0.01, and a marginally significant one for 0.01<p<0.10 
(Nemec & Nemec [T985> . As compared with other tests, the ran- 
domization test is more rigorous in rejecting peaks that might 
contain some oscillatory signal; on the other hand there is no 
doubt that the peaks that it finds significant represent an os- 
cillation (Muglach 2003 1. For the period P=27.4 days we got 
p=0.009 which implies an almost zero probability that the ob- 
served time series oscillations could have occurred by chance. 
From the extreme case of complete exclusion from the data 
analysis of the upper limits the result is that the period of 
P=27.4 days is still determined as a dominant period in a noiser 
periodogram with p=0.099. 

The gamma-ray emission is predicted to be produced via 
inverse Compton scattering of stellar photons by the relativis- 
tic electrons of the jet at each periastron passage (Taylor et al. 
1992 1996 1. Are the two gamma-ray peaks in the EGRET data 
indeed located around periastron passage ? To determine the 
orbital phase of the gamma-ray emission we have folded the 
data with the period P=26.496. As shown in Fig. |4] the two 
gamma-ray peaks occur at O r _ ra> ,=0.2-0.5. 

3. X-rays, and optical observations 

Variations in the structure of the accretion disk around a com- 
pact object in a X-ray binary system can be revealed by X- 
ray data when they show a transition of the spectral state from 
the high/soft state to the low/hard state (Tanaka 1997 1. In their 
high/soft state systems hosting black holes have a power-law 
component E~ r (Tanaka fl 9971 his Table 3) with T > 2. When 
the high/soft state evolves into a low/hard state the value of T 
for the power law becomes F ~ 1.6 (in this case for both X-ray 
binary systems containing black holes or neutron stars) (Tanaka 
119971 . The spectral transition from high/soft to low/hard corre- 
sponds to a change in the disk structure (Tanaka 1997 1. The 
radio emission from microquasars with steady jets is always 
correlated with low/hard states (Fender 2004 1. 
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Fig.l. By Tavani et al. 1998 EGRET light curve of 
2CG 135+01. 




Fig. 2. Result by the Phase Dispersion Minimization analysis. 
Note that the analysis provides the most probable period as 
a minimum (Stellingwerf 1978 1 and the vertical scale is here 
shown reversed. The dominant feature is at P=27.4+7.2 days. 



Changes of the hardness ratio in the source GRS 1915+105 
have shown to correspond to a drastic change of the inner ra- 
dius of the accretion disk (Belloni et al.]T997 1 followed by the 
onset of synchrotron radiation, first at infrared wavelengths and 
afterwards at radio wavelengths (Mirabel et al. 1998 1. 

Figure |5}top shows the resulting photon index F of RXTE 
data of LS I +61°303 published by Greiner and Rau ffflyft . 
As one can see, the range for T is 2.0-2.4 except at one point 
where the photon index is T=1.6. This point at T=1.6 is clearly 
below the average and therefore assumed to be affected by un- 
known systematic uncertainties. In the context of microquasars 
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Fig. 3. Gamma-ray observations of Fig. ^ folded with a pe- 
riod of 27.4 days. The Phase is computed with respect to to= 
2443366.775 JD (Gregorv l2"0"0"2l . The circles correspond to the 
data set "a", the full squares to the data set "b", the crosses to 
the data set "c" and the triangles to the upper limits. 
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Fig. 4. Gamma-ray observations of Fig. 1 folded with the radio 
period of 26.496 days. 



this value corresponds - like the change in the hardness ratio as 
discussed above - to a variation in the accretion disk followed 
by the ejection of matter. And indeed the simultaneous radio 
observations show (Fig. [3} Bottom) that the photon index value 
r=1.6 coincides with the onset of a radio outburst. 

The error bars in the photon index are large and there is only 
one critical point where the transition takes place. However, 
Taylor et al. (|l996l also noticed a hardening of the X-ray emis- 
sion at the onset of a new radio outburst. Moreover, Leahy et al. 



( 1997 ) report on two X-ray observations, both with power-law 
indices in the range 1.63-1.90 at orbital phases 0.2 and 0.42. 
While for the first value we can only speculate that it might 
be coincident with the first accretion rate peak (no simultane- 
ous gamma-ray observations are available), the second value is 
certainly coincident with the onset of a radio outburst (Fig. 3 in 
Harrison et al.|2000). 

In conclusion, the fit by Greiner and Rau (2001 1 suggests 
that 1) LS I +61°303 with his power-law component with F > 
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Julian Date - 2440000.0 

Fig. 5. By Greiner & Rau 2001 Top two panels: X-ray spectral 
fit parameters. Bottom two panels: X-rays and radio data. 



2 is mainly in the high/soft state typical for systems hosting 
black holes, and 2) it changes this state into a low/hard state just 
before/during (it is difficult to determine the exact time delay) 
the onset of a radio outburst. 

The suggestion - coming from the X-ray data - of a black 
hole as the accreting object is appealing. Only Punsley ( 1999 1 
discussed this possibility in detail and presented a model for 
the high-energy emission based on it. In the literature it has 
generally been assumed that LS I +61°303 hosts a neutron star 
because of an estimated mass of M v =1.5 M© by Hutchings & 
Crampton ( 1981 1. To derive that value Hutchings & Crampton 
( 1981 1 assumed a mass function /=0.02, an inclination i^70° 
(sin 3 i=0.8) and a mass for the companion star of M=10 Mo. 
Below, we will show that these values are only average values 
from a rather large possible range. 

Concerning the mass function / as noticed by Punsly 
( 1999 1 the observations of Hutchings & Crampton (1991) im- 
ply a value in the range 0.0028< / <0.043. 

Ultraviolet spectroscopy (Hutchings & Crampton (1981) 
indicates that the primary star is a main sequence B0-B0. 5 star 
(L~ 10 38 ) erg sec -1 with a possible range for the mass of 5-10 
MO. Values in the range 10-16.5 M© have been suggested by 
Punsly J1999I and values in the ranges 10-18 Mo by Marti & 
Paredes (fT995l >. 

Assuming 600 km s _1 as the critical rotational velocity for 
a normal B0 V star and that Be stars may rotate at a velocity not 
exceeding 0.9 of this value (Hutchings et al. 198 1 ), the very low 
limit for the inclination of the orbit compatible to the measured 
value Vsin i = 360 + 25 km s 1 (Hutchings & Crampton [l98U 
is of 38° (Massi et al. l2001i . On the other hand, the observed 
(Hutchings & Crampton 1 198 It shell absorption in the strong 



Balmer and He I lines for a disk sufficiently flat corresponds to 
a large inclination angle (i=90°) (Kogure 1969). 

In conclusion, the uncertainties of the parameters are 
0.0028< / <0.043 for the mass function, 38° <i< 90° for 
the inclination angle and 5 Mo < M< 18 M© for the mass. 
The chosen values of /=0.02, i^ 70° and M=10 M© are in- 
deed reasonable average values. However, already assuming 
i=38°, one obtains M A = 2.5 Mo. Assuming in addition M=18 
M© the result is M A =3.4 Mo, clearly a black hole. Finally, tak- 
ing /=0.043 we are faced with an even more massive object. 
We can conclude that on the basis of the present optical data 
the presence of a black hole in the system LS I +61°303 cannot 
be ruled out. 

4. Conclusions 

The emission of the gamma-ray source 2CG 135+01 is consis- 
tent with a periodic behaviour, with a period similar to the or- 
bital/radio period of the system LS I +61°303. The gamma-ray 
peaks, observed at two different epochs (Tavani et al. 1998}, re- 
main confined around the periastron passage, a fact which may 
imply that the most of the seed photons for Comptonization are 
indeed stellar photons (Taylor et al. U996> . 

The X-ray observations along an orbital cycle (Greiner & 
Rau 2001 ) show that the system remains always in a high/soft 
state, characterized by a photon index F > 2. Only at the onset 
of radio emission there happens a state transition, characterized 
by T ~ 1.6, to the low/hard state (Tanaka 1997, Fender et al. 
1998). These transitions between spectral states are related to 
the dramatic change in the structure of the accretion disk pre- 
ceeding the ejection of part of it into a jet (Belloni et al. 1 19971 
Mirabel etal. 1998 Fender etal. 1998). Good sampled simulta- 
neous X-ray and radio observations in different outbursts would 
be quite important to trace (by monitoring the photon index T ) 
the two accretion-ejection processes occurring in this peculiar 
source. 

Finally, we suggest that LS I +61°30 may host a black hole, 
because its high/soft state has a power law component with 
r >2, which is typical for systems hosting black holes. On the 
other hand the value of M A =1.5 M© - always quoted in the lit- 
erature - is related to average values for the inclination of the 
orbit, mass of the Be star and mass function. The ranges for 
these parameters presently available are so large that a massive 
black hole in LS I +61°303 cannot be excluded. A new deter- 
mination of all these values is desirable. 

Acknowledgements. I am grateful to Karl Menten, Jiirgen Neidhofer, 
Guidetta Torricelli, Leslie Hunt, Edward Polehampton and the un- 
known referee for their comments and suggestions. I acknowledge the 
data analysis facilities provided by the Staiiink Project which is run by 
CCLRC on behalf of PPARC and I am thankful to Jennifer Hatchell 
for the local help. 

References 

Belloni, T., Mendez,M., King, A.R., van der Klis, M. & van Paradijs, 

J. 1997, ApJ.,479, L145 
Fender, R., Corbel, S., Tzioumis, T. et al., 1999, ApJ, 519, L165 



Maria Massi: LS I +61° 303 in the context of microquasars 



5 



Fender, R. P 2004, in 'Compact Stellar X-Ray Sources', eds. W. H. G. 

Lewin and M. van der Klis, Cambridge University Press, in press 

|astro-ph/0303339| 
Gregory, P. C, & Taylor, A. R. 1978, Nature, 272, 704 
Gregory, PC, Peracaula, M, & Taylor, A.R. 1999, ApJ, 520, 376 
Gregory, PC. 2002, ApJ, 575, 427 
Gregory, PC, & Neish, C. 2002, ApJ, 580, 1 133 
Greiner, J., and Rau, A. 2001, A&A 375, 145 
Hall, T.A., Bond, I.H., Bradbury, S.M. et al. 2003, ApJ, 583, 853 
Harrison, F.A., Ray, PS., Leahy, D.A., Waltman, E.B., & Pooley, G.G 

2000, ApJ, 528, 454 

Hartman, R.C., Bertsch, D.L., Bloom, S.D., et al. 1999, ApJS, 123, 79 
Hutchings, J.B., Nemec, J.M., & Cassidy, J. 1979, PASP, 91, 313 
Hutchings, J. B., & Crampton, D. 1981, PASP, 93, 486 
Kaufman Bernado, M. M., Romero, G. E., & Mirabel, I. F. 2002, 

A&A, 385, L10 
Kogure, T. 1969, PASJ, 21, 71 

Leahy, D. A., Harrison, FA., & Yoshida, A. 1997 ApJ, 475, 823 
Leahy, D. A. 2001 A&A, 380, 516 
Maraschi, L., & Treves, A. 1981, MNRAS, 194, 18 
Marti', J., & Paredes, J. M. 1995, A&A, 298, 151 
Massi, M., Paredes, J.M., Estalella, R., & Felli, M. 1993, A&A, 269, 
249 

Massi, M., Ribo, M., Paredes, J. M., Peracaula, M., & Estalella, R. 

2001, A&A, 376, 217 

Massi, M., Ribo, M., Paredes, J. M., Garrington, S. T, Peracaula, M., 

& Marti, J. 2004, A&A, 414, LI 
Meier, D.L. 2001, ApJ, 548, L9 

Mirabel, I.F, Dhawan, V., Chaty, S., Rodriguez, L.F., Marti, J., 
Robinson, C.R., Swank, J., & Geballe, T.R. 1998, A&A, 330, L9 
Muglach, K. 2003, A&A, 401, 685 
Nemec, A. F, & Nemec, J. M. 1985, AJ, 90, 2317 
Paredes, J.M., Estalella, R., & Rius, A. 1990, A&A, 232, 377 
Paredes, J.M., Marti, J., Peracaula, M., & Ribo, M. 1997, A&A, 320, 
L25 

Paredes, J.M., Marti, J., Ribo, M., & Massi, M. 2000, Science, 288, 
2340 

Punsly, B. 1999, ApJ, 519, 336 
Stellingwerf, R.F 1978, ApJ, 224, 953 

Tanaka,Y. 1997, Accretion Disks-New Aspects, E. Meyer-Hofmeister 
and H. Spruits (Ed.), Springer- Verlag Berlin Heidelberg New 
York, 1. 

Tavani, M. 1995, in: The gamma-ray Sky with GRO and SIGMA, 

Dordrecht: Kluwer, 181 
Tavani, M., & Arons, J. 1997, ApJ, 477, 439 

Tavani, M., Kniffen, D., Mattox, J. R., Paredes, J. M., & Foster, R. S. 

1998, ApJ, 497, L81 
Taylor, A.R., & Gregory, PC 1982, ApJ, 255, 210 
Taylor, A.R., & Gregory, PC 1984, ApJ, 283, 273 
Taylor, A.R., Kenny, H.T., Spencer, R.E., & Tzioumis, A. 1992, ApJ, 

395, 268 

Taylor, A. R., Young, G, Peracaula, M., Kenny, H. T, & Gregory, 

P. C. 1996, A&A, 305, 817 
Zamanov, R.K. & Marti 2000, A&A, 358, L55 



